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Abstract 

We study the status of the CPT violating neutrino mass spectrum which has been proposed 
to simultaneously accommodate the oscillation data from LSND, KamLAND, atmospheric and 
solar neutrino experiments, as well as the non-observation of anti-neutrino disappearance in short- 
baseline reactor experiments. We perform a three-generation analysis of the global data with the 
aim of elucidating the viability of this solution. We find no compatibility between the results of 
the oscillation analysis of LSND and all-but-LSND data sets below 3a CL. Furthermore, the global 
data without LSND show no evidence for CPT violation: the best fit point of the all-but-LSND 
analysis occurs very close to a CPT conserving scenario. 
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I. INTRODUCTION 

The joint explanation of the oscillation signals observed in LSND in solar [3, IE El 
El IE E and atmospheric |E IE EE El | neutrino experiments, and in the KamLAND reactor 
experiment [I2j provides a big challenge to neutrino phenomenology. In Refs. 
was observed that the LSND signal could be accommodated with the solar and atmospheric 
neutrino anomalies without enlarging the neutrino sector if CPT was violated. Once such a 
drastic modification of standard physics is accepted, oscillations with four independent Am 2 
are possible, two in the neutrino and two in the anti-neutrino sector. The basic realization 
behind these proposals is that the oscillation interpretation of the solar results involves 
oscillations of electron neutrinos with Am^ < 10~ 4 eV 2 EE> while the LSND signal for 
short-baseline oscillations with Am 2 SND > 10 _1 eV 2 stems dominantly from anti-neutrinos 
(V^ — > V e ). If CPT was violated and neutrino and anti-neutrino mass spectra and mixing 
angles were different 

EE EE EE EE El 

both results could be made compatible in addition 
to the interpretation of the atmospheric neutrino data in terms of oscillations of both 
and z/ M with Am 2 tm ~ 10 -3 eV 2 0. 

In the original spectrum proposed, neutrinos had mass splittings Am 2 Q = Am^ <C 
Am^ = Am 2 tm to explain the solar and atmospheric observations, while for anti-neutrinos 
Am 2 tm = Amjj <C Am| x = AmL SND . Within this spectrum the mixing angles could be 
adjusted to obey the relevant constraints from laboratory experiments, mainly due to the 
non-observation of reactor !7 „ at short distances IE EE , and a reasonable description of the 
data could be achieved lg, |2l|, UM- m general, stronger constraints on the possibility of 
CPT violation arise, once a specific source of CPT violation which involves other sectors of 
the theory is invoked [23f . For a summary of recent theoretical work and experimental tests 
see, for example, Ref. [24] and references therein. 

On pure phenomenological grounds, the first test of this scenario came from the Kam- 
LAND 12| experiment since the suggested CPT-violating neutrino spectrum allowed to 
reconcile the solar, atmospheric and LSND anomalies, but, once the constraints from reac- 
tor experiments were imposed, no effect in KamLAND was predicted. The observation of 
a deficit in KamLAND at 3.5a CL clearly disfavoured these scenarios. Furthermore, Kam- 
LAND results demonstrate that V e oscillate with parameters consistent with the LMA v e 
oscillation solution of the solar anomaly This fact by itself can be used to set constraints on 
the possibility of CPT violation j^, 25, ~2^. Within the present KamLAND accuracy, how- 



ever, the bounds are not very strong because KamLAND data does not show a significant 
evidence of energy distortion. 

The present situation is that the results of solar experiments in v oscillations, together 
with the results from KamLAND and the bounds from other v reactor experiments show 
that both neutrinos and anti-neutrinos oscillate with Amg, Am^ < eV 2 . Adding 

this to the evidence of oscillations of both atmospheric neutrinos and anti-neutrinos with 
Am 2 tm ~ 1CT 3 eV 2 , leaves no room for oscillations with Am 2 SND ~ 1 eV 2 . The obvious 
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Figure 1: Post KamLAND CPT violating neutrino mass spectrum proposed in Ref. [~Lt| ]. 

conclusion then is that CPT violation can no-longer explain LSND and perfectly fit all other 
data 0. 

This conclusion relies strongly on the fact that atmospheric oscillations have been ob- 
served for both neutrinos and anti- neutrinos with the same Am^ tm . However, atmospheric 
neutrino experiments do not distinguish neutrinos from anti-neutrinos, and neutrinos con- 
tribute more than anti-neutrinos to the event rates by a factor ~ 4-2 (the factor decreases 
for higher energies). Based on this fact, in Ref. 0] & n alternative CPT- violating spectrum 
was proposed as shown in Fig. In this scheme only atmospheric neutrinos oscillate with 
Am atm and give most of the contribution to the observed zenith angular dependence of 
the deficit of /i-like events. Atmospheric dominantly oscillate with AmL SND which leads 
to an almost constant (energy and angular independent) suppression of the corresponding 
events. For low energies oscillations with Am^ eac can also be a source of zenith-angular 
dependence. The claim in Ref. [l^ was that altogether this suffices to give a good descrip- 
tion of the atmospheric data such that the scheme in Fig. ^ can still be a viable solution 



to all the neutrino puzzles. This conclusion was contradicted in Ref. |21[ by an analysis of 
atmospheric and K2K data. However, according to the authors in Ref. [17[ an important 
point to their conclusion was the consideration of the full 3z/ and 377 oscillations, while the 
analysis in Ref. 2l| was made on the basis of a 2z/+2l7 approximation. 

In this article we determine the status of the CPT violating scenario in Fig. ^ as an 
explanation to the existing neutrino anomalies. In order to do this, we perform a three- 
generation global analysis of the solar, atmospheric, reactor, and long-baseline (LBL) data, 
and compare the allowed parameter regions from this analysis to the ones required to explain 
the LSND data. We find that no consistency between the parameters determined by the 



1 This possibility was also discussed in the first version of Ref. |21j . 
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analyses of both data sets appears below 3a CL. 



II. NOTATION AND DATA INPUTS 



mf—rn/j and Am?. 



2 2 

mf — m A 



In what follows we label the states as in Fig. [T]with Am? , ,,, ; _,,,, ( 
We denote by U and U the corresponding neutrino and anti-neutrino mixing matrix |2j| 
which we chose to parametrize as 28] 



U 



• C13C12 S12C13 S13 ^ 

— S12C23 — S23S13C12 C23C12 — S23S13S12 S23C13 
\ S23S12 — Sl3C23 c 12 — S23C12 — S 13 S 12C23 C23 c 13 J 



(1) 



where c^- = cos 9ij and = sin fly and with an "over-bar" for the corresponding anti- 
neutrino mixing. In writing Eq. (0) we take into account that for our following description 
it is correct and sufficient to set all the CP phases to zero. 

In the anti-neutrino sector the reactor experiments provide information on 

the F e survival probability: 



1 — c^ 3 sin 2 26*12 sin 2 



A' 



m 



21 1 



- sin 2 26 



13 



-2 • 2 

c 12 sm 



1 - sin 2 2fl i3 sin 2 



AE 
Am^L 
AE 



i—2 -2 

+ s 12 sm 



4£ 



s 4 



+ C 



-4 

13 



sin 2 2fl 12 sin 2 



4£ 



Am\ 2 L 
4£ 



for Am\ x LjE < 1 
for Am 2 31 L/E > 1 



(2) 



In our analysis we include the results from the KamLAND 0], Bugey 20 1 and CHOOZ 
reactor experiments. For KamLAND we include information on the observed anti-neutrino 
spectrum which accounts for a total of 13 data points. Details of our calculations and sta- 
tistical treatment of KamLAND data can be found in Ref. 2^ . For reactor experiments 
performed at short-baselines we include the constraints from Bugey [2^ and CHOOZ 
which are the most relevant for Am 2 > 0.03 eV 2 and 0.03 eV 2 > Am 2 > 10~ 3 eV 2 , re- 
spectively. In our analysis of the CHOOZ data we include their energy binned data which 
corresponds to 14 data points (7-bin positron spectra from both reactors, Table 4 in Ref. [lii ]) 
with one constrained normalization parameter. For the analysis of Bugey data we use a total 
of 60 data points given in Fig. 17 of Ref. j^, where the ratio of the observed number of 
events to the one expected for no oscillations is shown for the three distances 15 m, 40 m, 
and 90 m. For technical details of our Bugey analysis see Ref. j^ . 

In the scheme under consideration the probability associated with the — > V e signal in 
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LSND is given by 

P LSND ee sin 2 2S LSND sin 2 (^ff^) = 4 ^ 2fl 13 sin 2 (^f^) , (3) 

where we have neglected terms proportional to Am^ which are irrelevant for LSND distances 
and energies. In Eq. (jSJ) we have introduced the notation 

Am LSND = A^3i , sin 2 26> LS nd = 4 sin2 2 #i3 , (4) 

which we will use in the presentation of our results. To include LSND we use the results 
of Ref. [3lj | . based only on the decay-at-rest anti- neutrino data sample, which has a high 
sensitivity to the oscillation signal. The y 2 is derived from a likelihood function obtained 
from an event-by-event analysis of the data |3J|. LSND has also studied the neutrino channel 
v n —> v e from decay-in-flight events. The full 1993-1998 data sample leads to an oscillation 
probability for neutrinos of (0.10 ± 0.16 ± 0.04)% [1], which, although consistent with the 
anti-neutrino signal, is also perfectly consistent with the absence of neutrino oscillations 
at LSND, as required in the CPT violating scenario. This fact is the first motivation and 
successful crucial test for the explanation of the LSND results by CPT violation. In view 
of the low statistical significance of the LSND neutrino signal we do not include it in the 
analysis. 2 

For the neutrino sector we use information from solar neutrino experiments and the 
K2K (33! LBL experiment. For the solar neutrino analysis we use 80 data points. We 
include the two measured radiochemical rates, from the chlorine Q and the gallium [5], 
0, 0] experiments, the 44 zenith-spectral energy bins of the electron neutrino scattering 
signal measured by the SK collaboration [3], and the 34 day- night spectral energy bins 
measured with the SNO detector. We take account of the BP00 j3 predicted fluxes 
and uncertainties for all solar neutrino sources except for the 8 B flux which we treat as a 
free parameter. For the relevant cases oscillations with Am^ are averaged out for solar 
neutrinos and the survival probability takes the form: 

Pl: M = 4 + 4 012) , (5) 

where P e 2 j /,sol (Am2 1 , #12) is the survival probability for 2v mixing obtained with the modified 
matter density N e — > c 2 3 iV e . 

The results of the analysis of the solar neutrino data 0] imply that Am^ has to be small 
enough to be irrelevant for the K2K baseline and energy, and the survival probability at 
K2K is 

pT = 1 - 4 (444 + 444) sin 2 • (6) 



2 We note, however, that because of a slightly different experimental configuration the data sample obtained 
from 1993-1995 had a higher sensitivity to the neutrino signal. From that data alone a 2.6a signal for 



v e oscillations was obtained [i^] which disfavoured the CPT interpretation. 
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In the analysis of K2K we include the data on the normalization and shape of the spectrum of 
single-ring /i-like events as a function of the reconstructed neutrino energy. The total sample 
corresponds to 29 events We bin the data in five 0.5 GeV bins with < E rcc < 2.5 
plus one bin containing all events above 2.5 GeV. The details of the analysis can be found 
in Ref. Q- 

Finally, the analysis of atmospheric neutrino data involves oscillations of both neutrinos 
and anti-neutrinos, and, in the framework of 3v+3V mixing, matter effects become relevant. 
We solve numerically the evolution equations for neutrinos and anti-neutrinos in order to 
obtain the corresponding oscillation probabilities for both e and fi flavours. In our calcula- 
tions, we use the PREM model of the Earth [36j matter density profile. We include in our 



analysis all the contained events from the 1489 SK data set [8j , as well as the upward-going 
neutrino- induced muon fluxes from both SK and the MACRO detector [10]. This amounts 
for a total of 65 data points. More technical descriptions of our simulation and statistical 
analysis can be found in Ref. j^. 

III. RESULTS 

Our basic approach to test the status of the scheme in Fig. ^ as a possible explanation of 
the LSND anomaly together with all other neutrino and anti-neutrino oscillation data is as 
follows. First, we perform a global analysis of all the relevant data, but leaving out LSND 
data. The goal of this analysis is to obtain the allowed ranges of parameters AmL SND and 
sin 2 26*lsnd as defined in Eq. (JU) from this all-but-LSND data set. We then compare these 
allowed regions to the corresponding allowed parameter region from LSND, and quantify at 
which CL both regions become compatible. 

In this approach we start by defining the most general x 2 f° r the all-but-LSND data set: 

Xall-but-LSND( Am L Am 31> #12, #23, #13^^21, A ^31, #12, #23, #is) = Xsol( Am 21, #12, #is) 
+ XK2K(A m 31, 6*23, #13) + XBugey+CHOOZ+KLAND^^l, Afn^ ,9 12 ,9 13 ) 

+ XATM( Am 21, Am L,#12,#23,#13|Am^, Am^,^ 12 ,^23,#13) • (7) 

Notice that in this comparison we have not included the constraints from the non-observation 
of 17^ — » I7 e transitions at KARMEN which, by themselves, disfavour part of the LSND 
allowed parameter region. The reason for this omission is that we want to test the status of 
the CPT interpretation of the 
LSND and KARMEN results 



SND signal using data independent of the "tension" between 



31 



We first focus on the parameters Am^ and #12 . These parameters are dominantly de- 
termined by solar neutrino data, which for any #13 prefer values of Am^ well below the 
sensitivity of atmospheric neutrino data. Therefore, solar data are mostly important to 
enforce the "decoupling" of the Am^ oscillations from the problem. In other words, the 
atmospheric neutrino analysis can be made without any loss of generality, in the standard 
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hierarchical approximation for neutrinos, neglecting the effect Am| but keeping the generic- 
3i> dependence on # 13 . Notice that, unlike in the CPT conserving case, in the relevant ranges 
of mass differences, #13 is not bounded by any "terrestrial" experiment. The dominant source 
of information on #13 is atmospheric data (and less important also solar data), and for this 
reason we consistently take into account this parameter in our analysis. Thus, after the 
marginalization over Am^ and 6*12 Eq. (jZJ) takes the form 

Xall-but-LSND( Am L #23, #13^7?^, Am^, 012 , 023, #ls) = Xsol.marg^M 

+ XK2K(Am3 X , 6> 23 , 613) + XBugey+CHOOZ+KLAND(Am 21 , Am| l3 9 12 , 0l 3 ) 

+ XATM(Am3i,#23,#i3|Am 21 , Am^, 0i 2 , 2 3, #13)- (8) 

Let us now discuss the information on Am 21 from reactor experiments. The observation of 
the V e deficit in KamLAND favours Am^ values near the best fit Am^ = 7 x 10~ 5 eV 2 . For 
such small values oscillations with Am^ have no effect for atmospheric neutrinos. Therefore, 
we will start by studying the case Am 2 21 < lfr 4 eV 2 in Sees. IfflAl and lllTBl We will relax this 
assumption in Sec. IIII C| where we investigate also a possible effect of larger values of Am^. 
Notice also, that the case of small Am^ is continuously connected to the CPT conserving 
scenario since it allows for CPT conservation in the "12" sector. For Am^ < 10 -4 eV 2 , one 
can easily marginalize over Am 21 and # 12 and Eq. (jHJ) further simplifies to 

XaU-but-LSND(Am3!, 6> 23 , 6*13 I Am^, 2 3, #13) = Xsol,margl2(#13) + xL-K^m^, 9 23 , 613) 

+ XBugcy+CHOOZ+KLAND,margT2( AM 31^13) + XaTm(A™31> #23, #13 I ^31 > #23, #13) • (9) 

Finally, we notice that for any value Am^ > 10~ 3 eV 2 the results from CHOOZ or, for 
larger values of Am 2 1; from Bugey, imply a strong limit on sin 2 2#i3 , and in order to obtain 
the V e disappearance observed in KamLAND 13 has to be small. Within this bound the 
results of the atmospheric neutrino analysis are almost independent of the exact value of 613 
and this parameter can be effectively set to zero in Xatm without any loss of generality. 

For the sake of concreteness we present the quantitative results corresponding to the 
normal ordering shown in Fig. ^ for both neutrinos and anti-neutrinos. We have verified 
that the conclusions hold also for the corresponding inverted orderings either for neutrinos 
and anti- neutrinos. Note that solar and atmospheric data require (Am^l <C (Am^j, and 
reactor and atmospheric data (and furthermore LSND) require |Am 21 | ^ lA^T 2 !]- For such 
hierarchies, the difference between normal and inverted schemes arises mainly from Earth 
matter effects in the propagation of atmospheric neutrinos and anti-neutrinos, and for the 
large values of | Am^ | required to explain the LSND signal Earth matter effects are irrelevant 
in the anti-neutrino channel. Within the present experimental accuracy, these effects are not 
important enough to lead to significant differences in the results of the atmospheric neutrino 
analysis for direct and inverted orderings (see for instance Ref. j^]). 
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A. Analysis of all-but-LSND data 



Using all the data described above except from the LSND experiment we find the following 
all-but-LSND best fit point: 
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(10) 



with Xaii-but-LSND min = 186.5 for 238 — 11 = 227 dof. 3 This can be directly compared to the 
corresponding analysis in the CPT conserving scenario: 



Am^ 
Am 2 ! 
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Am^ = 7.1 x 10~ 5 eV 2 
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(11) 



with Xaii-but-LSND min = 187 for 238 — 6 = 232 dof. We conclude that, allowing for different 
mass and mixing parameters for neutrinos and anti-neutrinos, all-but-LSND data choose a 
best fit point very close to CPT conservation and maximal 23 mixing. 

Next we illustrate the amount of CPT violation which is still viable. In order to do 
so we plot in Fig. El the allowed regions for the largest neutrino and anti-neutrino mass 
splittings Am 2 x versus Am 2 ! and the mixing angles #23 versus 823 and #13 versus #13 and (after 
marginalization with respect to all the undisplayed parameters). The different contours 
correspond to regions allowed at 90%, 95%, 99% and 3a CL for 2 dof (A^ 2 = 4.61, 5.99, 
9.21, 11.83, respectively). In general the regions are larger for anti-neutrino parameters as a 
consequence of their smaller contribution to the atmospheric event rates. In particular Am^ 
can take values below the region of sensitivity of CHOOZ. As a consequence the limit on 
#13 at high confidence level is very week. Our results for the allowed regions for the largest 
neutrino and anti-neutrino mass splittings show good qualitative agreement with the 2u 
analyses of Refs. (s[ . In particular we find that within the 2v oscillation approximation 
the atmospheric neutrino analysis rejects the CPT violating scenario at a level close to 4a 
(x 2 (Am 2 tm = Awlsnd) — Y 2 (Amj = Am 2 tm ) = 14) which is in reasonable agreement with 
the ~ 5cr rejection obtained in Ref. j2l|. As expected, the introduction of the 3u mixing and 
the reactor data leads to some quantitative differences in the size of the allowed regions. 



The 10 neutrino parameters shown in Eq. I|1U|) plus the free solar 8 B flux give a total of 11 fitted parameters. 
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Figure 2: Allowed regions for the largest neutrino and anti-neutrino mass splittings Am^ and 
Am|n and the mixing angles #23 and #23 , and #13 and #13 (after marginalization with respect to 
the undisplayed parameters) for Am^Ara^ < 10 -4 eV 2 (see text for details). The different 
contours correspond to the two-dimensional allowed regions at 90%, 95%, 99% and 3a CL from 
all-but-LSND data. The best fit point is marked with a star. 



From our results shown in Eqs. f|10j) . (jllj) and in Fig. |2]we conclude that current global 
neutrino oscillation data excluding LSND show no evidence for CPT violation, since the 
best fit point is very close to a CPT conserving scenario. However, from present data a 
sizable amount of CPT violation by neutrino parameters is allowed (for a recent discussion 
on the comparison with the limits existing on the K — K mass difference j3] see Ref. 
We note that it will be possible to significantly improve the limits on CPT violation in the 
neutrino sector by future experiments such as neutrino factories, see for example Ref. j^. 

Concerning LSND, we find that values of Am^ = Am 2 SND large enough to fit the LSND 
result do not appear as part of the 3a CL allowed region of the all-but-LSND analysis 
which is bounded to Am^ < 1.6 x 1CT 2 eV 2 . The upper bound on Am^ is determined 
by atmospheric neutrino data (and slightly strengthened by the reactor constraints). To 
illustrate the physics behind this result we show in Fig. |3] the zenith-angle distributions 
of various atmospheric data samples for "Point A" with the following parameter values: 
Am 2 x = 2.5 x 1CT 3 eV 2 , Am^ = 0.9 eV 2 , s 2 23 = s 2 23 = 0.5, s 2 13 = 0.05, sj 3 = 0.005, 
A774 < 10~ 4 eV 2 , and Am^ < 10~ 4 eV 2 . This point has been chosen to be compatible 
with the LSND result while keeping an optimized Xaii-but-LSND- As seen m ^ ne figure this point 
fails in reproducing the up-down asymmetry of muiti-GeV muons as a consequence of the 
angular-independence in the deficit of the anti-neutrino events. Furthermore, it predicts a 
too large deficit of up-going muon events near the horizon since I7 M oscillations with Am 2 x = 
0.9 eV 2 lead to the disappearance of z/ M 's even at those higher energies and shorter distances. 
For up-going muons the contribution from anti-neutrino events is only half of that from 
neutrino events. As a consequence this data sample is most sensitive to the anti-neutrino 
oscillation parameters. Both effects, the wash out of the up-down asymmetry and the deficit 
of horizontally arriving up-going muons, contribute in comparable amounts to the statistical 
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Figure 3: Zenith-angle distributions (normalized to the no-oscillation prediction) for the Super- 
Kamiokande e-like and /i-like contained events, for the Super-Kamiokande stopping and through- 
going muon events and for Macro up-going muons. The full line gives the distribution for the best 



fit of v u — > v T oscillations and CPT conservation: Am, 
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the expected distributions for typical LSND-compatible CPT violating cases with the following 
parameter values: Point A: Am^ = 2.5 x 1CT 3 eV 2 , Am^ = 0.9 eV 2 , s 2 ^ = s 2 3 = 0.5, sf 3 = 0.05, 
s 2 3 = 0.005, Am%x < 10~ 4 eV 2 , and Am^ < 10" 4 eV 2 ; Point B: Am^ = 2.5 x 10" 3 eV 2 , Afn| = 
0(eV 2 

and s\ 2 = 0.75. 
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0.5, s| 3 = 0.25, sf 3 = 0.05, s 2 13 = 0.005, Am^ < 10~ 4 eV 2 , Ara^ = 5 x 10" 4 eV 2 , 
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Figure 4: 90%, 95%, 99%, and 3a CL allowed regions (filled) in the (Am| = Am| SND , sin 2 26> L snd) 
plane required to explain the LSND signal together with the corresponding allowed regions from 
our global analysis of all-but-LSND data. The contour lines correspond to A% 2 = 13 and 16 (3.2a 
and 3.6a, respectively). 

disfavouring of the CPT violating scenario. 



B. Comparison of the all-but-LSND and the LSND data sets 

It is clear from these results that the CPT violation scenario cannot give a good descrip- 
tion of the LSND data and simultaneously fit all-but-LSND results. The quantification of this 
statement is displayed in Fig. |U where we show the allowed regions in the (Am^ = Am 2 SND , 
sin 2 26*lsnd) plane required to explain the LSND signal together with the corresponding al- 
lowed regions from our global analysis of all-but-LSND data. 

Fig. H] illustrates that below 3a CL there is no overlap between the allowed region of the 
LSND analysis and the all-but-LSND one, and that for this last one the region is restricted 
to Amjj = Am 2 SND < 0.02 eV 2 . At higher CL values of Am^ ~ 0(eV 2 ) become allowed - 
as determined mainly by the constraints from Bugey - and an agreement becomes possible. 
We find that in the neighbourhood of Am 2 x = Am 2 SND = 0.9 eV 2 and sin 2 2#lsnd = 0.01 the 
LSND and the all-but-LSND allowed regions start having some marginal agreement slightly 
above 3a CL (at A% 2 = 12.2). A less fine-tuned agreement appears at 3.3a CL (Ax 2 ~ 14) 
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for Am 2 31 = Am 2 SND > 0.5 eV 2 and sin 2 2fl LSND < 0.01 . 

Alternatively the quality of the joint description of LSND and all the other data can be 
evaluated by performing a global fit based on the total x 2 -f unc ti° n Xtot = Xaii-but-LSND + 
Xlsnd; an d a PPlyi n g a goodness-of-fit test. The best fit point of the global analysis is 
sin 2 2fl LSND = 6.3 x 10~ 3 and Am 2 31 = 0.89 eV 2 with xL,mm = 200.9. In the following we 
will use the so-called parameter goodness-of-fit 41, 4^|, which is particularly suitable to test 



the compatibility of independent data sets. Applying this method to the present case we 
consider the statistic 

-2 _ 2 _ 2 _ 2 

"X. Xtot,min Xall-but-LSND,min XLSND,min (12) 
= A Xall-but-LSND( b - f O + A XLSN D ( b - f > 



where b.f. denotes the global best fit point. The \ 2 °f Eq. (fl2|) has to be evaluated for 2 dof, 
corresponding to the two parameters sin 2 29i,p^p and Am 2 x = Am 2 SND coupling the two 
data sets: all-but-LSND and LSND (see Ref. [42] for details about the parameter goodness- 
of-fit). From Ax 2 u _ but _ LS ND = 12 - 7 and a Xlsnd = l - 7 we obtain x 2 = 14.4 leading to the 
marginal parameter goodness-of-fit of 7.5 x 10~ 4 . 



C. The effect of large Am^ 

Finally, we study whether the conclusions of the previous subsections could be affected 
by allowing for larger values of Am^, such that its effect can show up in the atmospheric 
neutrino data and improve the quality of the fit as suggested in Ref. (l7j |. In Fig. El we show 
the dependence on Am^ of the \ 2 obtained for the analysis of atmospheric and CHOOZ 
data, and for all-but-LSND data. In each curve we have marginalized with respect to the 
undisplayed variables subject to the condition Am 2 x > 1 eV 2 . For the sake of normalization 
we have subtracted in each case the corresponding Xmm cpt f° r the CPT conserving scenario. 

The figure shows that, indeed, considering only atmospheric+CHOOZ data there is an 
improvement (although mild) in the quality of the fit due to the effect of oscillations with 
larger values of Am 2 x . To illustrate this result we show in Fig.Olthe zenith-angle distributions 
of various atmospheric data samples for "Point B", which gives the lowest Xaii-but-LSND f° r 
a larger value of Am\ x : Am\ x = 2.5 x 10" 3 eV 2 , Am^ = C(eV 2 ), sj 3 = 0.5, s 2 23 = 0.25, 
sj 3 = 0.05, s 2 13 = 0.005, Am 2 x < 10" 4 eV 2 , Am 2 21 = 5 x 10" 4 eV 2 , and s\ 2 = 0.75. The figure 
shows that the main effect of Ami os cillations is to increase the number of contained e-like 
events, in particular sub-GeV 39, 4^, improving the fit for those events. However, the two 



main sources of discrepancy in the atmospheric fit in these scenarios - the small up-down 
asymmetry for multi-GeV muon-like events and the deficit of horizontally arriving up-going 
muons - remain a problem even when atmospheric anti-neutrino oscillations with the two 
relevant wavelengths are included. We conclude from this analysis that the claim in Ref. ^| 
of a possible improvement of the atmospheric neutrino fit due to the inclusion of the effect 
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Figure 5: A^ 2 = X^in c^T ~~ ^min cpt as a function of Am^j from the analysis of atmo- 
spheric+CHOOZ data (lower line) and from all-but-LSND (upper line) data (see text for details). 

of oscillations with larger values of Am^ is qualitatively correct for the contained events, 
although quantitatively relevant only for the sub-GeV e-like events. Moreover, we find that 
quantitatively the improvement in the fit is not enough to make the scenario viable. This 
conclusion is partially based on the bad description of the upward-going muon events in the 
CPT violating scenario, a fact which was overlooked in Ref. ["T^ . 

In other words, our results show that atmospheric neutrino data are precise enough to 
be sensitive to the anti-neutrino oscillation parameters, and it cannot be well described 
by a combination of neutrino oscillations with Am^ ~ 3 x 10~ 3 eV 2 and anti- neutrino 
two-wavelength oscillations with Am^ ~ 0(eV 2 ) and Am^ ~ few 10~ 4 eV 2 . 

The net effect in the global all-but-LSND analysis is that the improvement in the atmo- 
spheric fit is not enough to make the scenarios viable because it does not fully overcome the 
preference of smaller Am^ in KamLAND (even within their present limited statistics) as il- 
lustrated in the all-but-LSND curve in Fig. 03 The local minimum at Am 2 x = 7.1 x 10~ 5 eV 2 
corresponds to a point in the vicinity of the point where the LSND and all-but-LSND re- 
gions in Fig. |3] first meet. From the curve in Fig. |S] we see that the improvement obtained 
by moving to the minimum at Am^ = 5 x 10~ 4 eV 2 is only 0.5 units in % 2 . We conclude 
that higher Am^ values do not significantly affect the overall status of the CPT violating 
scenario. 



IV. CONCLUSIONS 

We have explored the possibility of explaining all the existing neutrino anomalies without 
enlarging the neutrino sector but allowing for CPT violation as described by the scenario 
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in Fig. ^ In order to do so we have performed a compatibility test between the results of 
the oscillation analysis of the LSND on one side and all-but-LSND data on the other in 
the framework of 3^+317 oscillations. Our main results are shown in Fig. HJ We find that 
the allowed regions for both data sets have no overlap at 3a CL. Alternatively, using the 
so-called parameter goodness-of-fit our results imply that the probability for compatibility 
between both data sets within this scenario is only 7.5 x 1(C 4 . 

The information most relevant to our conclusion comes from the atmospheric neutrino 
events. Our results show that, within the constraints imposed by solar and LBL neutrino 
data, and reactor anti-neutrino experiments, atmospheric data are precise enough to be 
sensitive to anti-neutrino oscillation parameters and cannot be described with oscillations 
with the wavelengths required in the CPT violating scenario. 

Furthermore, the global oscillation data without LSND show no evidence for any CPT 
violation. An analysis of the all-but-LSND data set allowing for different mass and mixing 
parameters of neutrinos and anti-neutrinos gives a best fit point very close to perfect CPT 
conservation. 
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